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I. INTRODUCTION 



During the last decade observational cosmology has advanced to the point of placing 
a number of constraints on our theoretical understanding of the universe. Observations 
of Type la Supernovae (SNe la) have become a critical probe of the expansion history of 
the universe |l| and recent measurements of the cosmic microwave background indicate an 
accelerated expansion of the universe occurred in an early inflationary epoch [2]. Current 
efforts to understand these observations within a mathematically consistent framework of- 
ten start with the ten or eleven dimensional spacetime manifolds motivated by superstring 
or M theory [3] Q. These theories often require a compactification phase to achieve the 
four dimensional spacetime manifold associated with the current universe and a means for 
describing the six dimensional manifold in terms of current observational constraints, i.e. 
they are not obviously visible. The observations indicate that the expansion of the universe 
is accelerating, at critical density, has a nonzero cosmological constant and is flat. These ob- 
servations also estimate that some 70% of the universe is composed of dark energy, another 
25% dark matter and the rest baryonic matter. The dark energy is often characterized by 
a negative pressure equation of state p = wp where the equation of state parameter for the 
exotic matter, w, lies near -1. One candidate for the exotic matter is a type of dark energy 
associated with a cosmological constant arising from zero point vacuum energy fluctuations 
[H[ 0that give a constant value of w = —1 or quintessence arising from a time varying scalar 
field which gives a dark energy that decreases with time — 1 < w < — 1/3 where dark energy 
density decreases with scale factor a(t) as p cx a^~ 3 ^ 1+w ^ or phantom energy corresponding 
to cases of exotic matter with w < — 1 Q Q . 

An understanding of dark energy can be searched for in string/M-theory and the result- 
ing generalizations of low energy gravity that are motivated by such inve stig ations. It is 
possible to utilize unusual cou pling s to matter fields | or a dilaton field |l0j ]. or coupling 



constants extra dimensions [I2I] or Kaluza Klein compactifications 3> or a time varying 
cosmological constant or dimensional reductions 0] or corrections to the pure gravity 
sector of the theory through interesting functions of the curvature 0] or Gauss-Bonnet 
motivated couplings |l8|or couplings to branes |l9|. The use of phantom fluids where the 
pressure and energy density are related by an equation of state given by p < — p; which for a 
scalar field can be achieved by using a negative kinetic energy term in the Lagrangian den- 
sity, is another promising avenue of development j2{]|2lj]. One demands that such theories 
give reasonable descriptions of the large scale four dimensional universe described by the 
Friedmann-Roberston- Walker (FRW) equations arrived at from the standard Einstein field 
equations. An insightful way to approach the dark energy problem is to utilize modifications 
of the Einstein field equations; to search for ways to accommodate the dark energy observa- 
tions within the context of higher dimensional or hig her order modifications to the theory 
or to use noncommutive geom etries and Casimir [22J energy or by equivalent scalar-tensor 
theories or moduli fields [23j. Nojiri and Carrol examine positive and negative powers of 
the scalar curvature in order to describe a stable inflationary phase and an accelerating 
late phase in geometric terms 24[ 25[ 26[. Within the framework of General Relativity the 
cosmic acceleration can be accommodated for with a phantom equation of state. Lovelock 
modifications consist of dimensionally extended Euler densities which include higher order 
curvature invariants in such a way that they lead to second order gravitational field equa- 
tions but are always total derivatives in four dimensions j2?| . The quadratic curvature terms 
are referred to as Gauss-Bonnet terms which come about in the low energy limit of heterotic 
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string theory 00. 

Here we consider a dynamic compactification of a (4 + d) -dimensional manifold to a maxi- 
mally symmetric manifold of dimension d and an expanding FRW space of dimension 4. We 
modify the gravity Lagrangain to first order in string tension a by a Gauss-Bonnet term that 
is quadratic in the curvature |3p|. Here we have considered these terms explicitly coupled 
to matter fields obeying some equation of state with the strength of each term treated as a 
parameter to be fixed by observations. It is known that when the internal space is a non- 
singular time independent compact manifold without boundary then there is no accelerated 
expansion unless the strong energy condition is violated and if the higher dimensional stress 
energy tensor obeys the strong energy condition than so does the lower dimensional stress 
energy tensor j3j| • Here we will consider the time dependent solutions of arbitrary dimension 
on a maximally symmetric manifold thereby avoiding the time independent nature of the 
compact space. In this paper we will investigate the Einstein-Gauss-Bonnet field equations 
on a manifold M n that undergoes dynamical compactification to M n = (R 1 x M 3 ) x S d . 
We will include matter by specifying an equation of state and applying the constraint equa- 
tions on the energy momentum tensor d^T^ = 0. Within this model we will extract the 
Friedmann-Robertson- Walker equations with the Gauss-Bonnet terms as correction terms 
to the standard 4-d equations. These equations are then solved pertubatively and analyzed 
within the context of the observational constraints and an equation of state. To understand 
the impact of higher order curvature terms on the observed cosmic evolution we will solve 
the Gauss-Bonnet corrected FRW equations that undergo a dynamical d-dimensional com- 
pactification to a maximally symmetric space. This paper is organized as follows. In Section 
[n]we present a general Einstein plus Gauss-Bonnet action in <i-dimensions where the field 
equations are calculated and the correction to the FRW equation is given. In section ITTT1 we 
prepare a perturbative solution for special cases including with and without the cosmologi- 
cal constant and calculate the corrected values of the Hubble constant and the deceleration 
parameter. In Section HVl we summarize our results. 



II. FRAMEWORK FOR FIELD EQUATIONS 

In this paper we will follow the notation of Paul and Mukherjee k| and Mohmmedi 
to express the Einstein-Hilbert action with an extra Gauss-Bonnet term as 



S = J d D x^[R-2X-eG] (2.1) 



where Q is a Gauss-Bonnet term (Q = R A bcdR ABCD — ^R A bR AB + R 2 )- A variation of the 
action (|2.1|) with respect to g A s produces the equation; 

G A C + \g A c + eg A c = ^T A C (2.2) 

where T is the 4 + d dimensional coupling constant and the Einstein and Gauss-Bonnet 
tensors, respectively, are 

G A C = R A C — \g A °R (2.3) 

Q A C = - (RbdefR BDEF — 4RbdR BD + R 2 ) 8 a c 

— (2RbdeaR BDEC + 2RR A C — AR D B R BA DC — AR B C R A B ) ■ (2.4) 
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Furthermore we will assume that the stress-energy tensor will be that of a perfect fluid, 
thus it is of the form 



T/ = diag[p(t),p(t),p(t),p(t),p d (t),...p d (t)\ (2.5) 

where p d (t) is the pressure of the higher dimensional compact manifold. 
We will choose a metric ansatz: 

dr 2 



ds 2 = -dt 2 + a 2 it) 



1 - Kr< 



+ r 2 (d6 2 + sin 2 6d(j) 2 ) 



+ b 2 (t) lmn (y)dy m dy n (2.6) 



where the extra dimension is defined to be maximally symmetric such that the Riemann 
tensor for 7 has the form R a bcd = ^{laclbd — ladlbcj ■ In agreement with current observations 
we will consider the usual 3 spatial dimensions to be flat (K = 0) and also demand that 
the extra dimensions be flat as well (k = 0) in agreement with the finding that one finds 
unphysical properties for p and p if k 7^ [14j . 

The metric leads to Riemann tensors of the following form (where both the dimensions 
are flat) 

RaOaO — a( i j Rabab = ^d 2 , -R m 0m0 = bb , Rmama = ddbb , R a bab = bv (2-7) 

where a, b are indices which run from 1...3 and m,n are indices which are in the extra 
dimensions, (obviously the symmetries as well R a bcd = Rbadc = —Rabdc = —Rbacd) The Ricci 
Tensor and Ricci Scalar are 

R 00 = 2,-J r A 1 R aa = 2d 2 + ad + d^- , R mm = 3— + (d- l)b 2 + bb (2.8) 
a b b a 

ii h hb n 2 \? 

R = Q- + 2d- + Qd-j + 6- + d(d -I)- (2.9) 
a b ab a z b z 

where d is the number of extra dimensions. 

The full Einstein tensor may be expressed as having an Einstein term and a Gauss- 
Bonnet term . We make the assumption that the extra dimensions compactify as the 3 
spatial dimensions expand as Q| 

where n > in order to insure that the scale factor of the compact manifold is both 
dynamical and compactifies as a function of time. With this ansatz, the non-zero elements 
of the Einstein tensor takes the form 

1 a 1 
G ° = - \6 + dn(dn-n~6)} — 

2 a 2 

= ^- 2 (2.H) 

cr 

G a a = (dn-2)--\[2 + dn{dn + n - 2)] ^ 



a 2 a 2 



■ i 

a a 



a a 2 



V2- + V3- (2.12) 

G m m = (dn-n-3)--i[6 + n(d-l)(dn-4)]4 

a 2 a z 



• i 

a a 



r}A-+r] b — . (2.13) 



a 'a 2 
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Note that when d = 4 the equations become the well known FRW equations in four dimen- 
sions. In 4 + d dimensions, the Gauss-Bonnet terms become 



r 



Qa a 



dn 



(d-l)n 



-(d-2)n[(d-3)n- 12] + 18 



12 



1 , 

-dn 

2 



(2.14) 



(d +l)(d- l)(d - 2)n 3 + 4(2 - d(d + l))n 2 - 4(d - 3)n + 



-2dn 



(d - l)n [(d - 2)n - 6] + 6 



aa 



d 4 



£2 -7 + £3 



•■ • 2 

aa 



(2-15) 



: (d- l)dn 2 



(d - 2)n [(d - 3)n - 8] + 12 



2(d- l)n 



(d - 2)n [(d - 3)n - 9] + U 



- 12 



aa 



•■ • 2 

aa 



a Liu, . , 

= e 4 - + e 5 — 2.16 

a 4 a 13 

where the constants rji and & depend upon the values of n and d as defined above. Note 
that if we allow d — > then the Gauss-Bonnet terms vanish for (? ° an d {? a a as one would 
expect in four dimensions. 



III. DYNAMIC COMPACTIFICATION 



For the case of both spaces having flat curvature, the D- dimensional Friedmann- 
Robert son- Walker (FRW) equations (eqn !2.2|) take the form 



_P_ 
2T 

p_ 

2T 

Pd_ 
2T 



as 



a a 

V2- + V3 — 

a a A 



Vi 

+ x 



+ e6 



d 4 



6-7 +e 



a 1 



■■ ■ 1 

aa 

a 6 



a a 
f)A- + ri b — + A 
a a z 



, _ d 4 dd 2 
+ e + 6 — 

a 4 a^ 



(3.1) 
(3.2) 
(3.3) 



As the scale factor of the compact manifold tends to zero with time, we will make the further 
assumption that the pressure in the extra d dimensions is zero (p^ = 0). Together with these 
Einstein equations, we also demand that the conservation equation hold for the stress-energy 
tensor (V^T^ = 0) or 



(a 3 p) + P ■ 



dV r ' "dt v ') b 
Using the assumption that b = 1/a™, this becomes 

£(«V)+4(« 3 ) 



+ d a 3 - (p + p d ) = 



(3.4) 



(3.5) 



which by simple algebra may be written in the more familiar form 



P 



3- (p + p) 
a 



0. 



(3.6) 



Note that this is simply a statement that dE = —P dV where we have defined an "effective" 
pressure (p) [3| which an observer constrained to exist only upon the "usual" 3 spatial 
dimensions would see as ^ 

P = p- -dn(p + p d ). (3.7) 

This effective pressure can be easily computed from the <i-dimensional FRW equations ()3.1|) - 
()3.3|) and is given by the relation 

(3.8) 



P 



1 



2T 3 Vl 



• 2 

a a 

2 - + - 

a a z 



+ 3 A 



■■ • 2 

aa 



7i ~T +72^" 



where we've defined the constants 



7i = 6 - T,dn £5 



72 = 6 - ^ dn (& + 



(3.9) 



Note that the above field equations (|3.1|) - (|3.3|) . (|3.8|) become the same as Mohammedi [3 
in the limit where e — > (ie. No Gauss-Bonnet terms). By redefining the coupling and 
cosmological constant in eqns. (J3.1|) and (J3.8)) . one recovers standard 4-D FRW cosmology 
for arbitrary values of n and d in the e — > limit. 



A. Asymptotic solution when A = 

We will first examine the case when the Gauss-Bonnet terms are small with respect to 
the Einstein terms so they may be treated as perturbations. We take the pressure in the 
extra dimensions, pa, to be zero so equation (|3.3|) becomes 

Pd d a 2 / a 4 ^ a d 2 \ , 

i£ = = 774- + 775-2 + e U4- + 6-- 3.10 
21 a cr \ a 4 or/ 

This equation may be solved perturbatively by assuming that a(t) ~ a (t) + e ai(t). Col- 
lecting terms in e, we get a zeroth order equation 

^ + ^4 = 0- (3-11) 

a 774 % 

This equation yields the solution 

a (t) = (/it + Z y)^/(r ? 4+r ?5 ) ( 312 ) 

where /x and 1/ are constants of integration and 774 + 775 7^ 0. 1 This solution can be used to 
find an expression for the zeroth-order energy density and effective pressure. We find 

Po_ _ 2 ViVl 1 (313) 

2T (7/4 + 775)2 

_Po_ 1 2 774 fa ~ 2775) 1 ,oi 4 N 

2T "2^ 1 (774 + 775)2 ' a ^ 4 +r, 5 )M 



1 ??4 + ?75 = demands a complex value of n which will not be considered here. 
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(3.15) 



Combining ()3.13|) and (|3.14p . we can eliminate ao and obtain an equation relating p and p 
of the form 

Po = - 3 1 - 2— p = wp - 
3 V W 

This equation gives a constraint on the possible values of n and d by choosing a value for u> 
and so it mimics an equation of state (EOS). At first glace it seems odd that the dynamical 
equations themselves impose a relation between p and p, but this happened because of our 
choice that the extra dimensions compactify as b(t) ~ l/a(t) n . Equation 13. 15l mav be solved 
for w as 

1 /3+(d- l)(dn-3)n 
w = - ' 



n - n . , 1 ' ( 3 - 16 ) 
n(d — 1) + 3 / 

Note, in the limit where n — > 0, w = 1/3 which is the relationship one would expect for 
a radiation dominated universe. Thus, the geometrical terms in the compactifacation are 
playing the same role as matter. Thus, by demanding that w have a physical value; one may 
use this relationship to restrict the choices of n and d. For instance if d = 7, then n must be 
less then 1/2 if w is demanded to have a physically reasonable value of between 1 and —2. 
The first order perturbation is of the form 



fli(t) + 



2 ai(t) + 

7/4 a 



3 ao +2 %o| 
a r]4 % 



ai(t) = -- (£4 ^ 4. 

T] 4 



dp Qp 
' 2 



Of) 



(3.17) 



Using the solution for ao as found above, the solution for a± can be found and is of the form 

Ql, horn (t) +a x>non (t) . (3.18) 

The homogeneous solution is given by 

ai, hom (t) = A(pt + ^WOh-hb) + B ^ t + ^-Wfa+w) (3 19 ) 

where A and B are constants of integration and a\ non is the non- homogeneous solution. The 
resulting non- homogeneous differential equation is of the Cauchy-Euler form (if you plug in 
for ao): 



^1, non (t) 



MP - 1; 



dlnonif) + 



p 2 (5((5 



^l,non(t) 



V4 



(pt + u) 1 ' nonv ' (pt + u) 2 
where we've defined the quantity 

^ (774 + 

The non-homogeneous solution is ai^ non (t) = C (pt + u) 13 ' 2 where 



k 



(pt + V 



4-/3 



c 



2r/ 4 



V5 



e 5 



(3.20) 



(3.21) 



(3.22) 



Thus the general solution for a(t) to first order is 



a(t) = (pt + vf 



1 + e 



.4 



B 



C 



(pt + v) (pt + u) 2 



(3.23) 
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We may use this form of a(t) to find the Hubble parameter and the acceleration for small e. 
We find 



H 

qH 2 



a 
a 
a 



Ho 



1 - 



P 



B 

(pt + u) 



+ 



1 - 



2e 

7T 



B 



2C 

(2(3 



+ 



3) C 



(pt + u) (13-1) fat + !/)'< 



(3.24) 
(3.25) 



where H and go Hq are the zeroth-order values of these parameters. Note that in the large 
time limit (t — ► oo) these terms tend to their zeroth-order values. Hence, we see that the 
Gauss-Bonnet contributions become vanishingly insignificant for late cosmological times as 
one would expect. 



B. Exact Solution when 775/774 = £4/^5 



We will now examine the equation for p^ (eqn 13.3)1 when the ratio of the coefficients of 
the Gauss-Bonnet terms is equal to the ratio of the coefficients of the Einstein terms or 



k 



(3.26) 



This equation holds for a wide variety of n and d. Notice that for the special case of 
parameters the equation for the <i-dimensional pressure (eqn I3.1()jl can be factored as 







a 774 a 2 



r/ 4 + e£ 5 ^ 



This equation can be solved exactly and has a solution 

a(t) = (fit + „)W(»J4-hb) 



(3.27) 



(3.28) 



This solution can be used to find an expression for the energy density and effective pressure 
from the remaining field equations. As was done above, we can eliminate a(t) and obtain 
an expression relating p and p. 



a 



'C(i-WC) 



(p - (p) 



2eap 



where we defined the constants 



a 



C 
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72 - 71- 



(3.29) 



(3.30) 



with w given by (|3.15jl . This equation relating p and p is inherently non-linear as one might 
expect from an exact solution to Gauss-Bonnet FRW cosmology. For late cosmological times, 
we expect a small Gauss-Bonnet contribution to the field equations (|2.2|) . In this small e 
regime, ()3.29|) reduces to 



(p - wp) - 



a 



2C(i-wo 



Q5-Cp) 2 ^0 



(3.31) 
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Notice that in the e — ^ limit (no Gauss-Bonnet terms), this expression yields a 4D equation 
of state (|3.15p . For a small but non-zero e, we have a small correction term to the FRW 
equation of state. This equation can be expressed as 

p = pw + ep 2 w' (3.32) 

where w' = (((8(w + 3)(d 2 + d - 3))/(d(d + 2))) for the special case of n — 2/d. 

We can also examine the regime where the Gauss-Bonnet term dominates over that of 
the Einstein term for early cosmological times. To examine this regime, we take the large e 
limit of ()3.29|) and obtain 

(p - (p) ~ — C(l - «VO(P - wp) - (3.33) 
to 

Notice that in the e — > oo limit, we obtain a linear relation between the density and pressure 
as in (j3.15j) but with the parameter w being replaced by (. 



C. Asymptotic solution when A = constant 



In the previous subsection, we examined Gauss-Bonnet FRW Cosmology for a vanishing 
cosmological constant. For a non-zero value of A, the higher dimensional pressure equation 
is of the form 



^ = 
2T 



V4- + ^ — + A + e + & — ? 



a a 



(3.34) 

a a* " ' 

Again, we may solve this equation perturbatively. The zeroth-order equation yields the 
solution 

a (t) = [A exp (6t)+B exp (-5 t)f (3.35) 

where 5 2 = A (rj^ + rj 5 )/rjf. As was done in the previous subsection, we may use the remaining 
zeroth-order field equations to obtain an expression relating p and p. We find the relation 



po = wp + 4/3(1 + 2w)\ 
The first-order perturbation is of the form 
V5 <V 



di(t) 



r/ 4 a 



ai(t) 



do 775 dn A 

3— + 2—4 + 4 — 

a r] 4 a z Q rj 4 



ai(t) 



1 

V4 



d dg 
' 2 

a a 



Using the solution for a (t), ()3.37|) becomes 



,V5 3 



P 4 



di(t) + ( 2^(3 ) di(t) + f3 2 ai (t) = -— 



V4 



A exp ((3 t) 



(3.36) 



a (3.37) 



(3.38) 



where we set B = as we are interested in an expanding universe and defined the constant 
P 2 = A/ (7/4 + 775). The solution for a\{t) can be found and is 



hom (t) + a linm (t) 
where the homogeneous and transient solutions are given by 

0>l,hcrm(f) = ao(t) 

a>i,ntm(t) = -X a>o(t) 



(3.39) 



(3.40) 
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where c\ and c 2 are arbitrary constants and 



7 2 = (Vlhl - 1) 
1 r/4/3 2 



X 



2775(^4 + ^5) 



^4 



(3.41) 



Notice that the sign of the quantity 7 2 determines whether the scale factor experiences 
oscillatory motion or exponential growth and decay. Explicitly, exponential growth and 
decay occur when 775/774 > 1 or 775/774 < — 1. In terms of n and d, when 



(d-l)(n-- 



3<^>0 



n(d-l)(2-dn) > 



(3.42) 
(3.43) 



Notice that ()3.42|) is trivially satisfied for all permissible values of n and d. However, ()3.43|) 
is only satisfied when n < 2/d; thus n = 2/d represents a critical point of the system for 
which value the scale factor behaves as if it is critically damped. For this critical point, we 
find that the coefficients dramatically simplify to 



Vl=V3=Vi = V5 



2 + d 
d 



772 = 



As in subsection IIII Al we may use this value of a(t) to calculate the Hubble and acceler- 
ation parameters in the limit where the Gauss-Bonnet terms are small with respect to the 
Einstein terms: 



H = - ~ H 

a 



1 + e 



qH" 



qoH 2 Q 



(7 - l/ftc^ 1 ^ 1 - (7 + l/P)c 2 e-^ +1 '^ 
l-e{ (1 - [7 + (1 - 1//3)] 2 ) c ie ^-W- 

(1 - [7 - (1 - l/P)] 2 ) c 2 e-^W 



(3.44) 



(3.45) 



Note that these terms are either exponentially decaying or growing. 



IV. CONCLUSION 

To summarize, we have found the form of the conformal factor, a(t), within the context 
of an action that is the sum of Einstein and Gauss Bonnet terms in 4 + d dimensions with 
a compactifing <i-dimensional flat Kaluza-Klein space for a variety of cases. By assuming 
that the extra dimensions compactify as a(t) ~ b(t)~ n , we showed that restrictions on the 
possible form of the equation of state are imposed. Using pertubative methods to solve 
the field equations, a general form for a(t) to order e was found for both a zero and a 
non-zero cosmological constant . We used these values of a(t) to calculate values for both 
the Hubble and the acceleration parameters (H(t),q(t)) and found a first order correction 
due to the Gauss Bonnet terms which can play the role of additional source terms. These 
solutions are consistent with having a dark energy era for reasonable values of n and d. An 
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interesting subcase is when 775/774 = £4/^5, for these values of n and d we are able to solve the 
d dimensional pressure equation exactly. This is consistent with and generalizes the results 
of both Paul and Mukherieejlil] and Mohammedi (TiJ . A special case happens if n = 2/d, 
when this is true then the equations simplify considerably as detailed above and may be 
indicative of an underlying symmetry. 

The Gauss Bonnet terms can only contribute significantly during the early phase of the 
universe as one would surmise from dimensional arguments since the Einstein terms go 
as R ~ t -1 / 2 but Q ~ R 2 ~ t^ 1 etc. It is apparent that the Gauss-Bonnet term does 
not contribute significantly to any Big Rip scenario and that the current EoS with exotic 
matter dominating over baryonic matter still gives rise to this singularity in the future. A 
scenario where a scalar field coupled to a Gauss-Bonnet term that avoids a Big Rip has been 
investigated by Nojiri, Odinstov and Sasaki j33|. However in the early universe the phantom 
and/or quintessence energies need not dominate over the Gauss-Bonnet contributions and 
the EoS in this regime can still be radiation dominated. 

This paper leads to several questions. What explicit mechanism might drive cementifi- 
cation and on what timescale does this occur? First how does the coupling constant (T) 
run as the extra dimensions compactify? A second issue is how the addition of these Gauss 
Bonnet terms would change the semi-classical states of FRW solutions to the Wheeler-de 
Witt equation. The general relation of time scales is also important since it is not uncommon 
to have all of the critical compactification processes complete on the order of a Planck time. 
These issues will be reported on in future papers. 
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